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A homothetic, static, spherically symmetric solution to 
the massless Einstein-Klein- Gordon equations is described. 
There is a curvature singularity which is central, null, bi- 
furcate, massless and marginally trapped. The space-time is 
therefore extreme in the sense of lying at the threshold be- 
tween black holes and naked singularities, just avoiding both. 
A linear perturbation analysis reveals two types of dominant 
mode. One breaks the continuous self-similarity by periodic 
terms reminiscent of discrete self- similarity, with echoing pe- 
riod within a few percent of the value observed numerically 
in near-critical gravitational collapse. The other dominant 
mode explicitly produces a black hole, white hole, eternally 
naked singularity or regular dispersal, the latter indicating 
that the background is critical. The black hole is not static 
but has constant area, the corresponding mass being linear 
in the perturbation amplitudes, explicitly determining a unit 
critical exponent. It is argued that a central null singularity 
may be a feature of critical gravitational collapse. 
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I. INTRODUCTION 

The discovery by Choptuik [Q of critical phenomena 
in gravitational collapse is widely regarded as the most 
unexpected success of numerical relativity. The now fa- 
mous mass scaling law and echoing behavior had not 
been anticipated by existing black-hole theory. Despite a 
subsequent theoretical framework involving dynamical- 
systems and renormalization-group methods (|,^, this 
is still not fully understood; see reviews by Choptuik 
and Gundlach [^. In particular, the exactly criti- 
cal space-time suggested by such numerical simulations 
is not known analytically. Thus it is of some interest to 
study solutions with similar properties. 

The numerical results suggest a self-similar critical 
space-time. For the massless Einstein-Klein- Gordon 
case originally studied by Choptuik, there is a dis- 
crete self-similarity with a characteristic echoing pe- 
riod Other matter models yield continuous self- 
similarity, which can be described by a homothetic vec- 
tor. A class of homothetic, spherically symmetric so- 
lutions to the massless Einstein-Klein-Gordon equations 
was given by Roberts and has been studied by vari- 
ous authors ||^-[Tl[|. This article describes another such 
solution, which is static. This solution, the Roberts class 
and its time-inverse form a two-parameter family which 



was actually given by Brady |8j , though in that reference 
a parameter was fixed to reduce to the Roberts class. 

Since these solutions are continuously rather than dis- 
cretely self-similar, they do not include the critical space- 
time itself. However, one may hope to find a perturbed 
space-time close to the critical space-time. Indeed, Frolov 
|ll| has argued that linear perturbations of the critical 
Roberts solution break the self-similarity from continu- 
ous to discrete. This article therefore studies the static 
solution and its spherically symmetric linear perturba- 
tions. There are indeed perturbations which break the 
continuous self-similarity to a certain periodic structure, 
actually falling short of discrete self-similarity for reasons 
explained subsequently. The static symmetry greatly 
simplifies the analysis. In particular, this allows the first 
explicit construction of perturbed black holes with finite 
mass from such a critically trapped space-time. 

The article is organised as follows. Section II gives the 
solution and describes the global structure of the space- 
time. Section III studies linear perturbations, explicitly 
giving the general solution to the linearized field equa- 
tions, then imposing appropriate boundary conditions 
and finding the relevant modes. Section IV describes 
a type of dominant mode which produces echoing. Sec- 
tion V describes the other dominant mode, which pro- 
duces black holes for a certain quadrant in amplitude 
space. Section VI shows that the other quadrants pro- 
duce white holes, naked singularities and regular disper- 
sal respectively. Section VII concludes with some specu- 
lations concerning critical collapse. 



II. THE SOLUTION 

The field equations for Einstein gravity with a mass- 
less Klein-Gordon field fixing the relative coupling con- 
stant, are 



R = 2V0 (g) V0 
V^^ = 



(1) 

(2) 



where R is the Ricci tensor and V the covariant derivative 
of the metric g. A spherically symmetric metric may be 
written locally in terms of the line-element 



(3) 



where dOF' is the line-element of the unit sphere and (r, 7) 
are functions of the null coordinates x^. Then the area 
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of the spheres is Awr"^, so that r is the areal radius. The 
field equations in these coordinates are |]l2| 

d±d±r - 2d±-id±r = -r{d±(l)f (4) 

rd+d-r + d+rd-r + e^'< /2^0 (5) 

r2a+9_7 - d+rd-r - e^^ /2 = -r^a+^S.^ (6) 

rd+d-4> + d+rd-(j) + d-rd+4> = Q. (7) 

It will also be convenient to use coordinates (p, r) defined 
locally by 

= ±eP^^. (8) 

The solution is given simply by 

r 7 4>^T (9) 

which is the case a — (3 = Q oi Brady |^ . The variables 
{X,Y) of Choptuik evaluate as (0,1/2). The line- 
element can be written as 

ds^ = e^P{d^'^ + 2dp^ - 2dT'^). (10) 

This is explicitly conformal to a flat metric times a unit 
sphere, so that a Penrose diagram can be given as in 
Fig.|l|. The metric is evidently static, with static Killing 
vector d/dr. It is also homothetic, or continuously self- 
similar, with homothetic vector 3/ dp. This means that 

La/dp9 = 2.9 (11) 

where L denotes the Lie derivative and the factor of 2 is 
conventional. It follows that the null vectors d/d£^^ are 
also homothetic vectors, where 

^±\n{±x^)^T±p. (12) 

In this sense the solution is doubly self-similar. Hyper- 
surfaces orthogonal to a homothetic vector may be called 
homothetic hypersurfaces. 

There is a central singularity, meaning r — there, 
at x^x~ = 0. This is a bifurcate null singularity, where 
the causal nature is defined by the conformal metric. It 
is therefore also marginally trapped, since i9-|-r = on 
x^ = 0. It is a curvature singularity, since the Ricci 
scalar is found to be just —r~^. The other conformal 
boundary corresponds to infinite r and constitutes null 
infinity, though the space-time is not asymptotically flat. 
Correspondingly, the mass or energy p^ , p^ 

m = (r/2)(l-5-i(Vr,Vr)) (13) 

evaluates as to = r/4, becoming infinite at null infin- 
ity. The singularity is massless in the sense that to = 
there. This confirms the view that massless central singu- 
larities lie at the threshold between generic trapped and 
naked singularities, which occur for to > and to < 
respectively [Oj. The space-time is summarized in Figj^. 
The (p, t) or^ coordinates take all real values, whereas 
the x^ coordinates have the range a;"*" > 0, x^ < 0. 



The future and past branches of the singularity occur for 
x^ — and a;+ = respectively, where d/dr is taken to 
be future-pointing. 




FIG. 1. Penrose diagrams of the space-time. The thick 
hne is the central (r — 0) singularity and the other conformal 
boundary is null infinity (r = co) . (i) The dashed lines are the 
static hypersurfaces (r constant) and the solid lines are the 
orthogonal homothetic hypersurfaces (p constant), (ii) The 
dashed lines are of constant x~ or ^~ and the solid lines are of 
constant or 5"^, which are also homothetic hypersurfaces. 

The solution just avoids being a black hole (or white 
hole) in the sense that there are no trapped surfaces, yet 
the singularity is marginally trapped. Also, the singular- 
ity just avoids being naked in the sense of Penrose p4[ , 
that it should lie to the future of some point and the past 
of another point on some observer's worldline. Likewise, 
an observer could see arbitrarily high but not infinite 
curvature from the future branch of the singularity. In 
summary, this is a new type of extreme space-time, lying 
at the threshold between black holes and naked singular- 
ities. 



III. LINEAR PERTURBATIONS 

Henceforth, spherically symmetric linear perturbations 
of the space-time will be studied, initially following the 
analysis of the critical Roberts solution by Frolov ||l^,|ll| . 
The linear perturbations studied here may be decom- 
posed by Laplace transformations into modes 

r ^ eP{l + er{T)e-''P) (14) 
7 = ej{T)e-''P (15) 
(I) = T + e4>{T)e-''P (16) 

governed by the complex frequency k. Assuming this 
form of the perturbations fixes the perturbative gauge 
freedom, apart from the zero point of t. The field equa- 
tions are linearized by removing terms o(e) in the per- 
turbation parameter, < e <C 1. The four independent 
components of the field equations yield three independent 



linearized equations: 

j={l-k/2f?-f"/4: (17) 

4> = ~P' ~ f/fc (18) 

r" + k^r + 2kj + 24>' = (19) 
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where /' = df/dr. The first two are solved for (7, <?i) 
exphcitly if f is known, so that the system reduces to a 
single fourth-order equation for f: 

f"" - 2{k^ -k + 2)P" + k^{k^ - 2fc + 4) = 0. (20) 

This linear ordinary differential equation can be solved 
via the auxiliary equation 

uj^ - 2(fc2 - fc + 2)cj2 + k^[k^ - 2fc + 4) = 0. (21) 

This is a biquadratic equation, so that the roots come in 
pairs, uj = ±0;^, ±cjb where 



UJA = (fc2 -2A: + 4)1/2 
ujb = k. 



(22) 
(23) 



Then the general solution for each k is 

r = A+e'^-*^ + A^e-'^^'^ + B+e"«^ + B^e"'^^^ (24) 

for constants A± , B± , except in special cases where roots 
coincide, when similarly standard solutions can be given. 
The special cases are (i) k = 0, when ujb = 0, (ii) k = 2, 
when oja = ojb, and (iii) k ^ 1 ± \/3i, when uja — 0. 
Weak boundary conditions will be taken as 



lini r 7^ lim r ^ 7^ 0. 



(25) 



If these conditions were not met, the perturbed r would 
vanish at the original infinity or become infinite at the 
original centre, respectively, either of which would be 
wildly different from the background space-time and 
therefore inconsistent with the linear approximation. 
Slightly stronger boundary conditions additionally im- 
pose 



3 lim r 

x^ — >±oo 



3 lim r 

x±^0 



(26) 



with similar justification. By writing a mode explicitly 
as 

= (-a;+x-)i/2 + e^(^+)(i-fe+-)/2(_a.-)(i-fc— )/2 (27) 

the weak boundary conditions can be seen to imply 

5ft(fc±w)<l (28) 

with cases of equality examined individually in the fol- 
lowing. For the lob modes this means 



a < 1/2 



(29) 



where a = 3?fc, b = whereas for the loa modes a little 
calculation yields 



a = 1 



&2 > 3 



(30) 



with UJA = 0^ ^ h"^ — 3. These modes are indicated 
in Fig.| 



CO. 



1/2 1 



FIG. 2. Relevant modes in the complex plane oi k = a + ib. 
The u)B modes fill a region while the uja modes form two lines. 



IV. ECHOING 

A common argument is that the relevant modes are 
those with a > and the dominant modes those with 
maximal a [||-^,|l^,|ll| . As depicted in Fig.^, this oc- 
curs for the UJA modes. One might also expect the 
dominant mode among these to be that with minimal 
since higher-frequency modes interfere destructively. 
This mode is given by fc = 1 ± v^i, case (iii) above: 



eP -f eAcoB{Vip + c). 



(31) 



The general solution for this mode would replace A with 
A -)- Cr, but if C 7^ the boundary conditions (|25| ) are 
not satisfied. The perturbed space-time is static and the 
perturbation 5r is periodic, with echoing in p of period 



A = 27r/V3 = 3.62. 



(32) 



which is within about 5% of the value 3.44 obtained nu- 
merically by Choptuik This suggests that the per- 
turbed space-time may be close to the critical space-time. 
However, the Klein-Gordon perturbation 5(j) vanishes for 
this mode. 

Note that the perturbed space-time is not discretely 
self-similar; despite 5r being periodic in p, r contains an 
exponential term in p. The situation may be understood 
schematically in terms of a phase-space picture. The uni- 
versality of the critical exponent suggests that the criti- 
cal solution is an intermediate attractor of codimension 
one ^-P. In other words, there is a critical hypersur- 
face in phase space, separating black-hole solutions from 
dispersal solutions. The background solution lies in this 
hypersurface, by the same argument as for the critical 
Roberts solution: the Roberts class, cut and pasted to 
flat regions as described by Brady 1^, describe either 
black-hole formation or dispersal, except for the delimit- 
ing critical case. Linear perturbations are represented by 
vectors from the background solution. One seeks a per- 
turbation tangent to the critical hypersurface, pointing 
towards the critical space-time. If the critical hypersur- 
face is curved, the perturbed solution will generally be 
outside the hypersurface and so exhibit exponential as 
well as oscillatory behavior. 

For the modes with k = 1 + ib, b^ > one finds 
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r = eP + eA cos{bp + c) cos(/3r + d) (33) 

so that the perturbations are periodic in both p and t. 
The variable X plotted by Choptuik is, to 0(e), 

i 1^ = - ^ sin(6p + c) sin(/3T + d) (34) 

which is explicitly sinusoidal in both p and r. Apart 
from the well known echoing in r ~ p, roughly sinusoidal 
oscillations in p at constant r do seem to develop in the 
numerical simulations, e.g. figure 3 of Choptuik 

V. BLACK HOLES 

The common argument P-p|, p^ , pl| that the dominant 
mode will occur for maximal ^k, as this is largest as 
p — oo, implicitly assumes that the other coordinate r 
is finite, which actually applies only where the singularity 
bifurcates. For the rest of the singularity, the explicit 
expression (^) shows that the dominant mode at 2;^ = 
instead occurs for maximal 3?(fc ±0;). As above, this 
maximal value is unity. The corresponding modes consist 
of those discussed in the previous section and the ujb 
modes with = 1/2. Non-linear perturbations may 
then be dominated by either or both, in a way which is 
unclear from a linear perturbation analysis. Henceforth 
the new modes will be discussed separately. 

The ujb modes with k = 1/2 + ib, b ^ 0, do not satisfy 
the boundary conditions (Eq), since r does not have a 
limit as 0, x^ 7^ 0. Thus the dominant mode is 

just LUg = k = 1/2. Then 

= {~x+x-f'^ + eB+{x+Y'^ + eB^{-x-f'^. (35) 

The locus of the centre r = can be determined by noting 
that it is a hyperbola in = (±x^)^/^ coordinates: 

Q = z+z- + eB+z+ + eB^z- . (36) 

Likewise, the locus of the trapping horizon, if any, is 
determined by 

d±r = ±{±x^)-^/H{Tx^y^^ + eB±)/2 (37) 

where d± = d/dx^. Here a trapping horizon |p^ , |l5| is 
a hypersurface foliated by marginal surfaces, 9+r — or 
d-r = 0, which respectively characterize black or white 
holes. These loci can be seen to be straight lines of con- 
stant x~ or x'^ respectively, or non-existent, depending 
on the signs of B± . These lines correspond to the asymp- 
totes 

z=F = ~eB± (38) 

of the above hyperbola, which passes through z'^ ~ z^ — 
0. The part of the hyperbola and its asymptotes inside 



z+ > 0, 2;~ > can be mapped back to the x coordi- 
nates, the loci depending qualitatively only on the signs 
of B±. This, combined with the nature of the background 
conformal boundaries, determines the qualitative global 
structure of the perturbed space-time. 

For the quadrant Bj^ < 0, B^ > 0, assumed in the 
remainder of this section, there is a black hole, with a 
region of trapped surfaces following a trapping horizon 
and preceding a spatial central singularity, as depicted 
in FigJ|. A mode with 5 = 0, < a < 1/2 is also de- 
picted for comparison, the analysis being similar. In that 
case, the trapping horizon is spatial, appears at zero area 
and evolves towards infinite area at conformal infinity. 
However, for the dominant mode, a remarkable feature 
is that the trapping horizon is null, occuring at constant 
x~ , with constant areal radius 

r = -e^B+B^. (39) 

This seems to be an example of an isolated horizon |T^ ] 
which is not a Killing horizon; the perturbed space-time 
is not static. It is also clearly an event horizon, gener- 
alized from the usual definition for asymptotically flat 
space-times. 



(i) / \ (ii) 




FIG. 3. Penrose diagrams of black holes produced by the 
ljb modes with 6 = 0, for (i) < a < 1/2 and (ii) the domi- 
nant mode a = 1/2. The thick curve is the central singularity 
and the other curve is the trapping horizon, between which 
lies the trapped region. This occurs for the quadrant B+ < 0, 
B- > 0. Dashed lines indicate conformal boundaries of the 
background space-time which are either outside or inside the 
perturbed space-time. 

A curious feature is that the past branch of the back- 
ground singularity is regularized by this perturbation: 
the X* coordinates break down at a;"*" = 0, but it can be 
shown to have finite area, finite mass m = r/2 and van- 
ishing Ricci scalar. (Essentially, 7 ^ 00 there, so norms 
V/ • V/ = —2e~'^'^d+fd-f tend to vanish, the exponen- 
tial term suppressing negative powers). This implies that 
it is a degenerate trapping horizon . Extension in- 

side this past event horizon seems, however, beyond the 
validity of the linear approximation; it is meaningless to 
have a linear perturbation of a region which does not 
exist in the background space-time. 

Since the mass or energy ([l3| ) is to = r/2 on a trapping 
horizon, there is a mass scaling law 

TO = -e'^B+B^/2. (40) 
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Thus the mass is exphcitly hnear in the perturbation am- 
pUtudes |ei3±| of both ingoing and outgoing modes, so 
that the critical exponent is unity. The renormaHzation- 
group argument , that the critical exponent is 

the reciprocal of the maximal JR/c, would not give this 
value unless is replaced by 3?(fc + oj), in accordance 
with the argument at the beginning of this section. This 
reflects the double self-similarity of the solution: the rel- 
evant homothetic vector at — is ±d/d£,^, rather 
than 8/ dp. 







black hole 


dispersal 




B+ 


naked singularity 


white hole 



FIG. 5. Summary of the effects of the dominant ujb mode, 
k — 1/2, in terms of the amplitudes B±. The background 
space-time lies at the intersection. 



VII. REMARKS 



VI. DISPERSAL, NAKED SINGULARITIES AND 
WHITE HOLES 

The other quadrants of the dominant lob mode, k = 
1/2, may be similarly analysed. The quadrant _B+ > 0, 
B_ < is a time-reverse of the black-hole case, yield- 
ing a white hole. For the quadrant B+ < 0, i?_ < 0, 
there is an eternally naked central singularity, with no 
trapped region outside the singularity. For the quad- 
rant B+ > 0, B- > 0, both branches of the background 
singularity are regularized: the area vanishes only at 
x+ = x~ = 0, a;^ = otherwise having finite area, finite 
mass m — r/2 and vanishing Ricci scalar, by a similar 
argument as for the black-hole case. Extension beyond 
these horizons again seems beyond the validity of the 
linear approximation. As far as can be judged from the 
unextended region, there is no singularity or trapped re- 
gion, with the perturbed solution describing regular dis- 
persal. The dispersal and naked-singularity solutions are 
depicted in Fig.^, the corresponding white-hole diagram 
being a time-reverse of Fig.^ii). 




FIG. 4. Penrose diagrams of space-times produced by the 
dominant lub mode, k = 1/2, as for Fig.^ii). (i) The quad- 
rant 5+ < 0, B- < 0, producing an eternally naked central 
singularity, (ii) The quadrant B+ > 0, B- > 0, describing 
regular dispersal. 

This indicates that, in the phase-space picture, the 
background space-time lies not only in the critical hyper- 
surface separating black-hole formation from dispersal, 
but also in another hypersurface separating black-hole 
formation from naked singularities, indicating a fourth 
region presumably describing white-hole decay. The pic- 
ture would be analogous to Fig.|^, with the critical space- 
time presumably lying roughly in the direction analogous 
to a third axis A for the dominant loa mode. 



The solution studied here describes a new type of ex- 
treme space-time, just on the verge of containing a black 
hole, white hole or naked singularity, without actually 
containing any of them. Each can be explicitly con- 
structed by linear perturbations, as can regular dispersal 
solutions, so that the space-time is also critical; it lies 
at the threshold between black holes, white holes, naked 
singularities and dispersal. The symmetry and simplicity 
of the solution suggest it as a standard example of an ex- 
treme critical space-time. As such, it may be useful in the 
context of cosmic censorship [^,0 . It should be noted 
that none of the perturbations studied here indicates a 
violation of cosmic censorship; the naked-singularity per- 
turbations involve an eternally naked singularity, not one 
forming from regular initial conditions. 

Various properties of the solution or its perturbations 
are analogous to those of near-critical numerical solu- 
tions in gravitational collapse. In particular, the dom- 
inant black-hole perturbation produces a black hole of 
constant mass and area, explicitly related to the pertur- 
bation amplitudes by a scaling law. Other perturbations 
produce echoing reminiscent of the numerical solutions. 

It may be conjectured that the numerically suggested 
critical space-time has a similar global structure to that 
of the solution described here, though this contradicts 
the analysis of Gundlach [||. Of course, the global struc- 
ture of a near-critical numerical solution is quite differ- 
ent, asymptotically flat with an initially regular centre. 
The critical and near-critical space-times should agree 
only close to the future branch of the singularity. A null 
central singularity in spherical symmetry is marginally 
trapped and therefore lies at the threshold of black-hole 
or white-hole formation, for the future-null and past-null 
cases respectively. An observer close to a future-null sin- 
gularity could see arbitrarily high but not infinite curva- 
ture, consistently with the numerical results. Likewise, 
a future-null singularity is not naked, in the sense that 
there is no future-causal curve from it, but is just on the 
verge of nakedness; to coin a phrase, a teasing singularity. 

Moreover, Christodoulou [|l8| showed analytically that, 
for the massless Einstein-Klein- Gordon system, there are 
non-generic collapses to a central null singularity and 
non-generic black holes preceded by a central null singu- 
larity, as depicted in Fig.^. This suggests the following 
picture of critical collapse: as a family of black-hole solu- 
tions approaches criticality, the trapped region is zipped 
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up between the trapping horizon and the centre. The 
closed zip is both a centre r = and, in a conformal 
sense, a trapping horizon d+r — 0, and therefore is a 
central null singularity. 




FIG. 6. Penrose diagrams of (i) generic and (ii) non-generic 
black holes formed from initial data on an outgoing null 
hypersurface (dashed line), according to the analysis of 
Christodoulou. The initially regular centre becomes a cen- 
tral singularity (thick curve) with the trapping horizon (the 
other curve) emerging as the singularity becomes spatial. In 
the generic case (i) the centre is instantaneously null, while 
in the non-generic case (ii) it is null for a finite time, during 
which it coincides with the trapping horizon in a conformal 
sense. A further limit of (ii) occurs when the singularity re- 
mains null for all time and there is no trapped region. 
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